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Abstract. We looked at the orbital motions of test particles according to the External Field Effect (EFE) predicted by 
the Modified Newtonian Dynamics (MOND) in the Oort cloud which falls in the deep MONDian regime (r ~ 50 — 150 
kAU). Concerning the interpolating function fl{x), we extensively used the forms fli = 1/(1 +x),/X2 = +^^)'''^iM3/2 = 
x/{l + j:-^/^)^/^. We integrated both the MOND and the Newtonian equations of motion in Cartesian coordinates sharing the 
same initial conditions. We considered both ecliptic and nearly polar trajectories, all with high eccentricities (e > 0.1). In 
order to evaluate the characteristic MOND parameters /x^ and Lg entering the problem, we used two different values (V = 220 
km s^' and V = 254 km s^') of the circular speed of the solar system's motion through the Galaxy; V allows to evaluate 
the Milky Way's gravitational field at the Sun's location. It turns out that EFE induces strong distortions of the Newtonian 
ellipses, especially in the ecliptic plane yielding more involved paths which span less extended spatial regions. 
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1. INTRODUCTION 

Modified Newtonian Dynamics (MOND) ^ was proposed by Milgrom in 1983 to explain the discrepancy between 
the observed kinematics of the exterior parts of spiral galaxies and the predicted one on the basis of the Newtonian 
dynamics and the baryonic matter detected from the emitted electromagnetic radiation (visible stars and gas clouds) 
Uj 2] without invoking exotic forms of still undetected non-baryonic Dark Matter. MOND postulates that for systems 
experiencing total gravitational acceleration A < Aq, with 

Ao = (1.2±0.27) X lO^'^ms^^ (1) 



A Amond = r. (2) 

r 

More precisely, it holds 

. ^Newton A 

A = — -— ,x=— ; (3) 

1 for X ^ I, i.e. for large accelerations (weak MOND regime, strong Newtonian regime), while /i(x) ^ x 
yielding eq. Q for x ^ 1, i.e. for small accelerations (deep MOND regime). The most widely used forms for the 
interpolating function are lUll^ 

M^) = (4) 

A'^C-) = (T^- (5) 

Such forms, and also another one, as we will see later, can be reduced to the following high-acceleration limit (x 1) 

H^l-kQX-". (6) 

Indeed, eq. (|4|i corresponds to A;o = 1, a = 1, while eq. Q corresponds to kQ = 1/2 and a = 2. 

In the framework of MOND, the internal dynamics of a gravitating system s embedded in a larger one S is affected 
by the external background field E of S, even if it is constant and uniform, thus implying a violation of the Strong 
Equivalence Principle: it is the so-called External Field Effect (EFE). In the case of the solar system, E would be 
Acen ~ lO"'" m s"2 because of its motion through the Milky Way USB]- 



2. EFE IN THE PLANETARY REGIONS OF THE SOLAR SYSTEM 



Perhaps the absence of explicit expressions for the action of EFE in the weak MOND regime, i.e. in the planetary 
regions of the solar system, may have lead to a misunderstanding in some sectors of the scientific community. 
Indeed, a researcher active in MOND writes: "MOND breaks down the Strong Equivalence Principle. This means 
that the acceleration of solar system's bodies depends indeed on the background gravitational field and not only on 
the tidal field. As shown by Milgrom, even if the external field was constant (and the tidal force vanishes), the internal 
acceleration would depend on the external field. Claiming that Acen is irrelevant is only valid if the field equation were 
linear" He/she also adds that "for trans-Neptunian objects and planets, one can ignore the Acen " Another researcher 
working on MOND tries to go in deeper details by writing: "For the main planets, the acceleration is much larger than 
Ao (the order of magnitude of the EFE), and the effect is negligible [...] The EFE maintains a constant direction in the 
planet revolution, and its effect cancels out. " Such statements are Ukely expressions of a widely diffuse belief about 
EFE in solar system. 

In fact, it has been shown in that a constant and uniform acceleration E = ExX + Ey.y + E,z directed along a fixed 
direction n and with E ~ Acen ~ Aq does induce non-zero long-period, i.e. averaged over one orbital revolution, effects 
on the Keplerian orbital elements of a planet. Moreover, the resulting perihelion precessions of the inner planets would 
be 4 — 6 orders of magnitude larger than the present-day limits on the recently estimated non-standard perihelion rates 

a- 

Milgrom [ 10] recently made a step forward by explicitly working out EFE in the deep Newtonian regime. It turns 
out to be equivalent to the action of a distant, localized body X || 1 1] 

Ax^-Jfr + 3Jf{r-nx)nx, (7) 

where 

(8) 

'•x 

is the so-called tidal parameter of X and «x is a constant unit vector in the direction of X. The acceleration of eq. (|7J 
comes from the following quadrupolar potential lfT2|] 

Ux~—[r^~3ir-hxf]. (9) 
The result by Milgrom can be obtained with the replacement 



where —q is the MOND quadrupolar parameter. Latest data from Saturnian perihelion ifTslfTill yield ifTTll 

-q>0.2 (11) 

with a maximum value 

-^max = 0.4-1. (12) 

However, the hypothesis that a planetary-sized body X really exists in the distant regions of the solar system should 
be regarded as a serious competitor because its confrontation with the planetary perihelion rates 1 11] yields plausible 
values for the distance of X for various choices of its mass, in agreement with previous theoretical predictions for X 
made to explain certain features of the Kuiper Belt and the Oort cloud. 



3. ORBITS OF OORT COMETS IN MOND 

Moving to the deep MOND regime in the remote periphery of the solar system, let us define the following quantities 

Acen ^ , /I T\ 

^ = ^>1, (13) 



X f dfi 

dx 



(14) 



lig = nir]), (15) 
L, = L(rj). (16) 



In the weak acceleration regime, for 

r»r,T7->/2, (17) 
i.e. in the Oort cloud ifTsll . the action of EFE is different, so that the total acceleration felt by an Oort comet is ll6l [To|] 



GM ( 2 2^ ^''^ 



z 

Note that, since the ecliptic longitude and latitude of the Galactic Center are about Agc ~ 180 deg, jScc ~ ^6 deg, 
EFE is directed along the X axis of the ICRF, i.e. the barycentric frame in which the motion of solar system's objects 
are usually studied. 
Concerning L, we have 

1 +x 1 +x 

X 1 

M2= , ^ L2 = T— ^, (20) 

X 1 

^^/2= (1+^3/2)2/3 ^ ^3/2 = Y^- (21) 

The form jU3/2 was proposed in fiot]. The value of Tj depends on the Galactic field at the solar system's location which 
can be obtained from its centrifugal acceleration 

''icen — (22) 

where Y is the speed of the Local Standard of Rest (LSR) and /? = 8.5 kpc is the Galactocentric distance. The standard 
lAU value for the speed is V = 220 km s^', but recent determinations fT^ obtained with the Very Long Baseline 
Array and the Japanese VLBI Exploration of Radio Astronomy project yield a higher value: V — 254 ± 16 km s^^ 
Thus, TJ ranges from 1.5 to 2.3. 



3.1. Ecliptic orbits 

3.1.1. Case \lT,j2 

We will, now, consider an Oort comet whose Newtonian orbit covers the entire extension of the Oort cloud. It has 
semimajor axis a = 100 kAU and eccentricity e = 0.5, so that its perihelion is at 50 kAU and its aphelion is at 150 
kAU; for the sake of simplicity, we will assume it lies in the ecliptic plane. Its Newtonian orbital period is fl, = 3 1 .6 
Myr We will, first, use tj = 2.0 {V = 254km s"') and ;U3/2, so that 

= 0.82 (23) 
Lg = 0.25 (24) 

Left panel of Figure [T] depicts the integrated Newtonian (dashed blue line) and MONDian (dash-dotted red line) orbits 
for the same initial conditions for —3Pb < / < 0. In the right panel of Figure [T] we show the case Tj = 1.5 {V — 
220 km s^^) yielding 

Hf, = 0.75 (25) 
Lg = 0.34. (26) 

The MOND trajectories are not closed and are much less spatially extended that the Newtonian ones; the overall 
shrinking of the orbit is more marked for the standard value of the LSR circular speed (right panel of Figure[T]i. 

Such an effect is particularly notable for highly elliptical Newtonian orbits, as shown by Figure |2] for e = 0.9 and 
—Pb < f < 0. Note also how wildly the MOND trajectory changes during one Keplerian orbital period. 
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FIGURE 1. Numerically integrated orbits of an ecliptic Oort comet with a = 100 kAU, e = 0.5, Pb = 31.6 Myr. Dashed blue 
line: Newton. Dash-dotted red line: MOND with 1X^/2- Left panel: Tj = 2.0 corresponding to V = 254 km s^' . Right panel: r] = 1.5 

corresponding to V = 220 km s^'. The initial conditions are xq = a{l — e),yo = zo = 0,xo = 0,yo = na^ l^'^o ~ 
span of the integration is —3f\, <t<0. 
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FIGURE 2. Numerically integrated orbits of an ecliptic Oort comet with a = 100 kAU, e = 0.9, = 3 1 .6 Myr. Dashed blue line: 
Newton. Dash-dotted red line: MOND with jl^/i- Upper panel: T] = 2.0 corresponding to K = 254 km s^'. Lower panel: T] = 1.5 

corresponding to V = 220 km s^'. The initial conditions are xq = a{l — e),.vo = = 0,io = 0,^0 = nfly^ T^'^" ~ ^' 
span of the integration is — Pj, < f < 0. 



3.1.2. Case ^2 



In Figure|3](e — 0.5) and Figure|4](e = 0.9) we use /X2. The values of the MONDian characteristic parameters are 



for T] = 2.0, and 



for Tj = 1.5. 



fi2 = TTTT. 1 = 2-0 (V = 254 km s"') 
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FIGURE 3. Numerically integrated orbits of an ecliptic Oort comet with a = 100 kAU, e = 0.5, P\, = 31.6 Myr. Dashed blue 
line: Newton. Dash-dotted red line: MOND with /i2- Left panel: 77 = 2.0 corresponding to V = 254 km s^'. Right panel: rj = 1.5 

corresponding to V = 220 km s^'. The initial conditions are xq = a{l —e),y\) = zo = 0,xq = 0,yo = nay^ j^'^o ~ 0- The time 
span of the integration is —3f\, <t<0. 



Also with such a form of the interpolating function fi, orbits that are highly eccentric in Newtonian dynamics are 
confined to much smaller spatial regions in MOND and experience high-frequency variations over one Keplerian 
orbital period (Figure |4|i. However, for jj.2 the largest extension of the MOND trajectory occurs for V — 220 km s^' 

(TJ = 1.5). 



3.1.3. Case /Xi 

Let us, now, examine the case /ij. In the left panel of Figure |5] we show the trajectory due to it of the Oort comet 
with e = 0.5 over 3Pb for T] — 2.0 which implies 

Hg = 0.67 (31) 
Lg = 0.33. (32) 



The case 77 — 1.5, yielding 



= 0.60 (33) 
Lg = 0.39, (34) 



is shown in the right panel of Figure|5]for —F\, < f < 0. The case of highly elliptic orbits (e — 0.9) is more intricate, as 
shown by Figure |6] The MOND paths resemble confuse clouds confined within a small spatial region. 
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FIGURE 4. Numerically integrated orbits of an ecliptic Oort comet with a = 100 kAU, e = 0.9, = 3 1 .6 Myr. Dashed blue line: 
Newton. Dash-dotted red line: MOND with \X2. Upper panel: r\ = 2.0 corresponding to V = 254 km s^'. Lower panel: f] = 1.5 

corresponding to K = 220 km s^'. The initial conditions are xq — a(\ — e),yo = zo = 0,xo = 0,3)0 = na^J~^^,ZQ = 0. The time 
span of the integration is —P\, < ? < 0. 




FIGURE 5. Numerically integrated orbits of an ecliptic Oort comet with a = 100 kAU, e = 0.5, P\, = 31.6 Myr. Dashed blue 
line: Newton. Dash-dotted red line: MOND with /ii. Left panel: 77 = 2.0 corresponding to V = 254 km s^'. Right panel: r\ = 1.5 

corresponding to V = 220 km s^'. The initial conditions are xq = a{l — e),.vo = zq = 0,xo = 0,yo = na^ I^jZo = 0. The time 
spans of the integration are — 3Pb < f < (left panel) and — < ? < (right panel). 



The general features common to all the pictures shown may have consequences on the interaction of the Oort-like 
objects lying close to the ecliptic with passing stars ifTsIl by reducing their perturbing effects and, thus, also altering 
the number of long-period comets launched into the inner regions of the solar system, the number of comets left in 
the cloud throughout its history. Indeed, in the standard picture, the comets moving along very (Newtonian) elongated 
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FIGURE 6. Numerically integrated orbits of an ecliptic Oort comet with a = 100 kAU, e — 0.9, = 3 1 .6 Myr. Dashed blue line: 
Newton. Dash-dotted red line: MOND with /ij. Upper panel: r) = 2.0 corresponding to V = 254 km s^'. Lower panel: r] = 1.5 

corresponding to V = 220 km s^'. The initial conditions are xq = a(l — e),yo = zo = 0,io = 0,3)0 = nay^ j^'^O ~ 
span of the integration is —0.9f\, <t<0. 



orbits may come relatively close to a star of mass suffering a change in velocity Ay which approximately is OlSll 

Av= — , (35) 

where is the star's velocity with respect to the Sun and d is the distance of closest approach with the Oort object. 
Moreover, less elongated orbits would also reduce the perturbing effects of the Galactic tides. 



3.2. Nearly polar orbits 

Let us, now consider the case of orbits showing high inclinations / to the ecliptic. For space reasons we will only 
show some cases. In Figure|7]we depict the sections in the coordinate planes of an orbit with a = 66.6 kAU, e — 0.92, 
7 = 81 deg for jUi and 77 = 1.5. 

The case of H2 and Tj = 1.5 is illustrated in Figure[8] 

It turns out that also in this case the MOND orbits are not closed, but the shrinking is now less marked than in the 
ecliptic plane. 



4. CONCLUSIONS 

The structure and the dynamical history of the Oort cloud, subject to EFE in deep MONDian regime, may be altered 
with respect to the standard Newtonian picture because highly eccentric orbits are not allowed in the ecliptic plane by 
MOND which, on the contrary, tends to strongly shrink them. As a consequence, one may speculate that the number 
of long-period comets launched in the inner parts of the solar system should be reduced because of the less effective 
perturbing actions of nearby passing stars, interstellar clouds and Galactic tides. Out of the ecliptic the situation is 



different because, although distorted with respect to the Newtonian case, the MOND orbits tend to occupy larger 
spatial regions than in the ecliptic. 
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FIGURE 7. Sections in tiie coordinate planes of the numerically integrated orbits of an Oort comet with a = 66.6 kAU, e = 0.92, 
7 = 81 deg. Dashed blue line: Newton. Dash-dotted red line: MOND with iL\,r\ = 1.5 (V^ = 220 km s^'). The initial conditions are 
XQ = 40 kAU,yo = 30 kAU,zo = 5 kAU,io = -23 kAU Myr-\jo = -15 kAU Myr-i,zo = -15 kAU Myr^i. The time span of 
the integration is — Pj, < f < 0. 
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FIGURE 8. Sections in tiie coordinate planes of the numerically integrated orbits of an Oort comet with a = 66.6 kAU, e = 0.92, 
7 = 81 deg. Dashed blue line: Newton. Dash-dotted red line: MOND with /i2> i] = 1.5 (V^ = 220 km s^'). The initial conditions are 
XQ = 40 kAU,yo = 30 kAU,zo = 5 kAU,io = -23 kAU Myr-\jo = -15 kAU Myr-i,zo = -15 kAU Myr^i. The time span of 
the integration is — Pj, < f < 0. 



